Effects of Interactions on the Critical Temperature of a Trapped Bose Gas 
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We perform high-precision measurements of the condensation temperature of a harmonically-trapped atomic 
Bose gas with widely-tuneable interactions. For weak interactions we observe a negative shift of the critical tem- 
perature in excellent agreement with mean-field theory. However for sufficiently strong interactions we clearly 
observe an additional positive shift, characteristic of beyond-mean-field critical correlations. We also discuss 
non-equilibrium effects on the apparent critical temperature for both very weak and very strong interactions. 
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The effect of inter-particle interactions on the Bose- 
Einstein condensation temperature of a dilute gas has been a 
topic of theoretical debate for more than fifty years, since the 
pioneering work of Lee and Yang fl]|2). In a uniform system 
there is no interaction shift of the critical temperature T c at 
the level of mean-field theory. However, consideration of the 
correlations between particles which develop near the critical 
point leads to the conclusion that repulsive interactions en- 
hance condensation, i.e. shift the condensation temperature 
above the ideal gas value T c ° OiJHt 

Ultracold atomic gases offer an excellent testbed for fun- 
damental theories of Bose-Einstein condensation and related 
many -body phenomena 1T31 Il4l . However, in these systems 
the problem of the interaction shift of T c is even more com- 
plex because they are produced in harmonic traps. In this case, 
at least for weak interactions, the T c shift is dominated by an 
opposing mean-field effect, which reduces the critical temper- 
ature lfl5l . Within experimental precision, previous measure- 
ments fl6l - [T8ll were consistent with the mean-field theory and 
could not discern the effects of critical correlations. 

In this Letter, we report on high-precision measurements of 
the T c shift in a potassium ( 39 K) gas with tuneable interactions 
Ifl9ll20l , We employ a Feshbach resonance ED to extend the 
previously explored range of interaction strengths and elim- 
inate several key sources of statistical and systematic errors. 
This allows us to clearly reveal the long-sought beyond-mean- 
field effects on the critical temperature. We also examine 
the stringent requirements for equilibrium T c measurements, 
which are violated in the regimes of either very weak or very 
strong interactions. In non-equilibrium gases we observe evi- 
dence for 'super-heated' condensates which survive at an ap- 
parent temperature above the equilibrium T c , suggesting that 
strong dissipation can stabilise the coherent condensed state. 

Historically, most theoretical work focused on a spatially 
uniform gas, and for several decades there was no consensus 
on the functional form, or even on the sign of the T c shift (for 
an overview see e.g. 171 [T0HT2ll ). It is now generally believed 
that the shift is positive and to leading order given by J71|8]: 

AT r i /q a 

-^«1.3anV3«1.8-, (1) 

where AT C = T c — X 1 ", a > is the s-wave scattering length, 
n the particle density, and A the thermal wavelength at tem- 



perature T®. The positive AT C implies that condensation oc- 
curs at a phase space density below the ideal gas critical value 
of nX 3 = C(3/2) ~ 2.612 (where £ is the Riemann function). 

For a harmonically trapped gas, the condensation tempera- 
ture is defined for a given atom number N, rather than for a 
given density n. For an ideal gas, k^,T® — HQ [N/^S)] 1 ^ 3 , 
where Q is the geometric mean of the trapping frequencies 
along the three spatial dimensions, and £(3) « 1.202. This 
corresponds to a phase space density in the centre of the trap 
equal to the uniform system critical value, n(0)A 3 = f (3/2). 
The interaction shift of the critical point can be expressed ei- 
ther as AT C (N) (for comparison with theoretical literature) or 
as AN C (T) (for easier visualisation, as in Figs. [TJ and |2|. 
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FIG. 1: (Color online) Opposing effects of interactions on the critical 
point of a Bose gas in a harmonic potential V(r). Compared to an 
ideal gas (dotted blue line) with the same T c , repulsive interactions 
reduce the critical density, but also broaden the density distribution 
(solid red line). Mean-field theory (dashed line) captures only the 
latter effect, and predicts an increase of the critical atom number N c 
at fixed temperature T, equivalent to a decrease of T c at fixed N . 

The two opposing effects of repulsive interactions on the 
critical point of a trapped gas are illustrated in Fig. [T] where 
we sketch the density distribution at the condensation point 
for an ideal (dotted blue line) and an interacting (solid red 
line) gas at the same temperature. In the spirit of the local 
density approximation, the critical density should be reduced 
by repulsive interactions. However, interactions also broaden 
the density distribution. For weak interactions the latter ef- 
fect is dominant, making the overall interaction shift AN C (T) 
positive, or equivalently AT C (N) negative. 

The negative T c shift due to the broadening of the density 
distribution in a harmonically trapped gas can be calculated 
using mean-field (MF) theory, which neglects the reduction of 
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the critical phase space density implied by Eq. ([T]i (see dashed 
line in Fig.[TJ). This approach gives ifTSl : 



AT C 
TP 



-3.426 



Ao 



(2) 



The dominance of the negative MF shift of T c over the positive 
beyond-MF one goes beyond the difference in numerical pre- 
factors in Eqs. Q and ([TJ. At the condensation point, in a non- 
uniform system only the central region of the cloud is close to 
criticality, which reduces the net effect of critical correlations 
so that they are expected to affect T c only at a higher order 
in a/\ . The MF result of Eq. (|2| should therefore be exact 
at first order in a/ An. There have been several attempts to 
theoretically combine the effects of MF repulsion and beyond- 
MF correlations on T c for a harmonically trapped gas E2U261 . 
but no consensus has been reached beyond the expectation that 
the additional beyond-MF shift should be positive. 

Previous measurements |fT6Ul8l . consistent with Eq. (|2| 
within experimental errors, were performed for a/ An rang- 
ing from 0.007 [18| to 0.024 [17|. We explore the range 
0.001 < a/Ao < 0.06, using the 402.5 G Feshbach resonance 
in the \F, mp) = |1, 1) lower hyperfine state of 39 K 11271 . As 
in ||20ll28l , we produce 39 K condensates in a crossed optical 
dipole trap which provides a close to isotropic trapping poten- 
tial. Near the bottom of the trap, ui/2tt = 75 — 85 Hz for our 
measurements with atom numbers A ps (2 — 8) x 10 5 . 

To measure the critical point we prepare a partially con- 
densed cloud, fix the optical trap depth, and let the atom num- 
ber decay towards A c through inelastic processes; meanwhile, 
elastic collisions redistribute particles between condensed and 
thermal components, and the temperature remains essentially 
constant ll28l . Initial preparation of clouds with various con- 
densed fractions is done at a = 135 ao, where ao is the Bohr 
radius. We then adjust a to the desired value by ramping the 
Feshbach field, and wait for an a-dependent hold time ihoid 
before releasing the gas from the trap and measuring its mo- 
mentum distribution through absorption imaging after 19 ms 
of time-of-flight (TOF). In the last part of the paper we discuss 
the strict requirements on the relationship between ihoid, the 
elastic scattering rate 7 e i, and the relevant atom-number decay 
time t for the measurements of T c to faithfully reflect equilib- 
rium properties of the gas. For now we focus on presenting 
the measurements which we trust to be in equilibrium. 

In addition to extending the a/Ao range, the Feshbach reso- 
nance provides us with two experimental advantages essential 
for both the precision and accuracy of our AT C measurements: 

(1) For each measurement series at a given a and Ao, we 
concurrently take a reference measurement with a different a, 
same ui and very similar A, hence very similar Ao. Specifi- 
ically, for the reference point we choose a small a such that 
a/Ao ~ 0.005. We thus directly access the small T c shift 
due to the difference in a/Ao, and essentially eliminate all 
a-independent systematic errors that usually affect absolute 
measurements of T C (N, Q, a). These include uncertainties in 
the absolute calibration of A and Q, as well as the additional 
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FIG. 2: (Color online) Determination of the critical point and dif- 
ferential interaction shift, (a) An absorption image of a cloud with 
0.14% condensed fraction, and the corresponding (azimuthally av- 
eraged) cut through the column density. The gas was prepared at a 
large scattering length, a — 274 ao, but the interactions were turned 
off in TOF. (b) Condensed (No) versus thermal (N 1 ) atom number 
for two concurrently taken series with a — 56 ao (blue circles) and 
a — 274 ao (black squares). Note that all points correspond to con- 
densed fractions below 2 %. The data is scaled to the same temper- 
ature (T = 240 nK) and shows the shift of the critical point in the 
form AN C (T). The solid point corresponds to the image shown in 
(a). Solid lines show the extrapolation to iVo = 0, necessary to accu- 
rately determine N c . (c) N' is plotted versus iV 2/5 for the same data 
as in (b), showing more clearly the extrapolation procedure. 



T c shifts due to finite-size effects ifTJI and the small anhar- 
monicity of the trapping potential 1 20 1 . 

(2) We home in on the critical point by turning off the inter- 
actions during TOF. To do this we quickly (in < 2 ms) ramp 
the Feshbach field to 350 G immediately after the release of 
the gas from the trap. This minimises the expansion of small 
condensates and allows us to reliably detect condensed frac- 
tions as small as ~ 10~ 3 (see Fig. [2}. 

Fig. [2] illustrates our differential measurement. Here a = 
27 '4 ao, Ao « 10 4 ao, and a = 56 ao for the reference series. If 
the two series had identical A c values, we could directly read 
off the differential AT C (N). To correct for the small (few %) 
difference in A c we apply the ideal gas scaling, T c oc A 1 / 3 , 
to the reference series. The second-order error in AT C due to 
the small (< 2 %) T c shift at a/A ~ 0.005 is much smaller 
than our statistical error bars. For visual clarity, in Fig.[2]we 
instead scale to the same temperature and display AN C (T). 

In Fig. [2jb-c) we show the relationship between the con- 
densed (Nq) and thermal (N 1 ) atom number near the critical 
point [29 1. The rise of Ao in Fig. [2jb) is not simply vertical 
because the thermal component in a partially condensed gas 
is not saturated at A c 1281 ; one can also see that this effect is 
more pronounced at higher a. It is therefore essential to care- 
fully extrapolate A' to the A = limit in order to accurately 
determine N c . We extrapolate using N' = N c + SoNq', 
with the non-saturation slope So(T,ui,a) calculated with no 
free parameters following [28 1. 

In Fig. [3] we summarize our equilibrium measurements of 
the interaction shift AT C /T®. We took data with a range 
of atom numbers, Aw (2 — 8) x 10 5 (corresponding to 
TJ? « 180 — 330 nK), in order to verify that our results de- 
pend only on the interaction parameter a/Ao. The dashed 
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blue line shows the MF result of Eq. which agrees very 
well with the data for a/Ao < 0.01. For larger a/Ao we ob- 
serve a clear deviation from this prediction. All data points 
are fitted well by a second-order polynomial (solid red line), 
AT C /T C ° = 61 (a/A ) + & 2 (a/Ao) 2 , with 61 = -3.5 ±0.3 and 
62 = 46 ± 5 [30 1 . From a theoretical point of view this func- 
tional form might not be exact, for example one might also 
expect small logarithmic corrections (see e.g. |24|). However 
within our error bars such corrections are not discernible. 

The value of b\ is in excellent agreement with the MF pre- 
diction of —3.426 lfl5l . The value of 02 strongly excludes 
zero, and its sign is consistent with the expected effect of 
beyond-MF critical correlations. These measurements pro- 
vide the first clear observation of beyond-MF effects on the 
transition temperature of a harmonically trapped gas. 
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FIG. 3: (Color online) Interaction shift of the critical temperature. 
Data points were taken with TV w 2 x 10 5 (blue circles), 4 x 10 5 
(black squares), and 8 x 10 (red triangles) atoms. The dashed line is 
the mean- field result AT C /T C ° = -3.426 a/X . The solid line shows 
a second-order polynomial fit to the data (see text). Vertical error 
bars show standard statistical errors. Horizontal error bars reflect the 
0.1 G uncertainty in the position of the Feshbach resonance. 

To conclude this part of the paper, we assess the systematic 
errors in our measurements. In general interactions increase 
the kinetic energy of thermal atoms during TOF; this results 
in an a-dependent error in T which does not cancel out in our 
differential measurements. This error is minimized by fitting 
the high-energy wings of the thermal distribution (excluding 
the central thermal radius from the fit) ||3T| . We also turn the 
interactions off at the beginning of TOF, but the reduction of a 
is gradual over « 2 ms. We measured the difference between 
(apparent) T with interactions 'on' and 'off during TOF to be 
approximately linear in a/Ao, an d about 4% for a = 400 ao 
and Ao ~ 10 4 ao. By varying the time at which we turn off a, 
we estimate our residual error to be 1 — 2 % at a/Ao = 0.04. 
This estimate is supported by numerical simulations which 
reproduce the experiments quite well. Additionally, interac- 
tions modify the initial in-trap momentum distribution. This 
reduces the apparent T because the positive chemical poten- 
tial preferentially enhances population of low-energy states. 
We numerically estimate this effect to be also approximately 
linear in a/Ao, an d about —2 % at a/Ao = 0.04. Fortuitously, 



the in-trap and in-TOF effects partially cancel, resulting in a 
net error in AT C /T C ° of at most ±1 % at a/A = 0.04. 

In the rest of the paper we discuss the equilibrium condi- 
tions required for our measurements, and the non-equilibrium 
effects revealed when they are violated. 

In general, a system with continuous dissipation can only 
be 'close to' thermodynamic equilibrium. For an atomic gas, 
the proximity to equilibrium depends on the dimensionless pa- 
rameter 7 c iT, which measures the relative rates of elastic and 
inelastic processes. In practice the j e i required for equilib- 
rium measurements also depends on the measurement preci- 
sion. We measure N c to about 1 %, so we require that the gas 
continuously (re-)equilibrates on a timescale r corresponding 
to only 1 % atom-loss. We thus require about 100 times higher 
7ci than one would naively conclude by taking the 1 /e lifetime 
of the cloud as the relevant timescale. Equilibration is usually 
considered to take about 3 collisions per particle [32]; for all 
measurements shown in Fig. [3] we made sure that -f c \T > 5. 
All our data also satisfy the condition <hoid > t > 1 /ui, nec- 
essary for global equilibrium to be established. 

An interesting question in its own right is what happens if 
we violate these stringent equilibrium criteria. In Fig. |4ja) we 
show measurements with N w 4 x 10 5 atoms (Ao ~ 10 4 ao), 
extending beyond the equilibrium region shown in Fig. [3] 
We still show only measurements that satisfy thoid > t > 
1/uj and 7ei£hoid > 5, so that there is nominally enough 
time for global equilibrium to be established. However if 
7ciT is not large enough, the elastic collisions cannot 'keep 
up' with the continuously present dissipation. The result- 
ing non-equilibrium effects can thus not be eliminated by 
simply extending <hoid> but are an intrinsic property of the 
system. In Fig. |4jb) we plot 7 c it, based on calculated 7 c i 
11331 and t measured near the critical point. Individually, 
7d ~ 0.7 — 1000 s~ 4 and t 2ms — Is vary vastly as a 
function of a {^ c \ increasing and t decreasing), but the break- 
down of equilibrium appears to occur at very similar values of 
7eir in the low- and high-a limit. 

Non-equilibrium phenomena necessarily depend on addi- 
tional factors such as the initial conditions, and we therefore 
do not expect our quantitative results to be universal. The de- 
tails of non-equilibrium dynamics will be a subject of future 
experiments, and here we discuss only qualitative trends. 

In the small-a limit we observe a smooth rapid rise of the 
apparent T c above the equilibrium curve (and hence above T c ° 
for a — > 0). We can qualitatively understand this effect within 
a simple picture. In this regime, losses are most likely dom- 
inated by one-body processes which equally affect Nq and 
N' . The net effect of equilibrating elastic collisions would 
therefore be to transfer atoms from the condensate to the ther- 
mal cloud. However the dissipation rate is too high compared 
to 7ei, and so No remains non-zero even after the total atom 
number drops below the equilibrium critical value iV c (i.e. the 
measured T c is above the equilibrium value). Note that strictly 
speaking T is not defined out of equilibrium, but the absolute 
value of the observed effect is sufficiently small that an equi- 
librium distribution function fits the data very well and pro- 
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FIG. 4: (Color online) Non-equilibrium effects, (a) AT C /T C ° for 
N ~ 4 x 10 5 atoms is determined following the procedure which 
assumes equilibrium (as in Fig. [2}. At both very low and very high a 
the apparent T c deviates from the equilibrium curve, (b) Equilibrium 
criteria (see text): 7 c ir (solid squares) is the number of elastic col- 
lisions per particle during 1 % atom-loss; J e i/Q = 1 (open circles) 
marks the onset of the hydrodynamic regime. 



vides a good measure of the energy content of the cloud. 

Our measurements in the large-a limit suggest that the ini- 
tial breakdown of equilibrium again results in condensates 
surviving above the equilibrium T c . However the physics in 
this regime is much richer, with several potentially competing 
effects requiring further investigation. For example, three- 
body decay affects Nq and N' differently, the thermal com- 
ponent is far from saturation ||28l , and the gas also enters the 
hydrodynamic regime, 7 e i/a) > 1 [see Fig. Qb)]. 

In conclusion, we have performed high-precision studies 
of the effects of interactions on Bose-Einstein condensation 
of a trapped atomic gas. In the regime where equilibrium 
measurements are possible, our most important observation 
is the clear deviation from mean-field behaviour for suffi- 
ciently strong interactions. The additional positive shift of the 
critical temperature is a clear signature of the condensation- 
enhancing effect of critical fluctuations. These measurements 
should provide motivation and guidance for further theoret- 
ical studies of this difficult problem. We also studied non- 
equilibrium condensation phenomena, for both very weak and 
very strong interactions. Further study of these effects should 
prove useful for understanding condensation in intrinsically 
out-of-equilibrium systems, such as polariton gases. 

We thank M. Holzmann and J. Dalibard for useful discus- 
sions, and M. Kohl for comments on the manuscript. This 
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